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NEW CLASSES OF NON-CONVOLUTION INTEGRAL 
EQUATIONS ARISING FROM LIE SYMMETRY 
ANALYSIS OF HYPERBOLIC PDFS 

MARK CRADDOCK AND SEMYON YAKUBOVICH 


Abstract. In this paper we consider some new classes of inte¬ 
gral equations that arise from Lie symmetry analysis. Specifically, 
we consider the task of obtaining solutions of a Cauchy prob¬ 
lem for some classes of second order hyperbolic partial differential 
equations. Our analysis leads to new integral equations of non¬ 
convolution type, which can be solved by classical methods. We 
derive solutions of these integral equations, which in turn lead to 
solutions of the associated Cauchy problems. 


1. Introduction 

The theory of continuous symmetry groups of systems of partial dif¬ 
ferential equations was developed by Lie in the last decades of the 
nineteenth century. A symmetry is a transformation which maps so¬ 
lutions to other solutions. Lie developed a method for systematically 
computing all continuous symmetries of a given system of differential 
equations. Excellent modern accounts may be found in the books by 
Olver [IB] and those of Bluman and his coauthors, such as the volume 

i. 

Symmetries are powerful tools and they allow us to solve a wide 
variety of problems. For example, we may compute fundamental so¬ 
lutions of parabolic equations from a knowledge of their symmetries. 
See the papers Ei.iia. 0. [10] for details, together with many exam¬ 
ples and applications. Boundary value problems may also be solved 
by the method of group invariant solutions. The aforementioned book 
[3] contains an extensive discussion of this topic. Symmetries are also 
essential to the study of conservation laws. A chapter of [IB] is devoted 
to this. 

A new way of applying Lie symmetries was introduced in [7] and de¬ 
veloped extensively in [B] . The idea is to construct an integral operator 
that maps test functions to solutions of the PDF by integration against 
a symmetry solution. The purpose of this paper is to show how the 
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method can lead to new classes of non-convolution integral equations, 
which may be solved by a combination of classical techniques. 

We hrst provide an example to illustrate how the method works. 
We then turn to a class of second order linear hyperbolic equations 
and derive some new integral equations which arise in the solution of 
Cauchy problems associated to these equations. We then solve these 
equations using a novel combination of integral transform methods. 

The outline of the paper is as follows. In Section 2 we introduce the 
method that forms the basis of our analysis. In Section 3 we determine 
the symmetries of a class of hyperbolic PDEs of the form Uu = u^x + 
f{x)u. We show that there are nontrivial symmetries when / satishes 
one of three families of Riccati equations. In Section 4 we study the 
equation Uu = y'^Uyy+yUy+y'^u ( which is equivalent to Uu = Uxx+e^^u, 
with y = e^) and show how the symmetries lead to a solution of the 
Cauchy problem for this equation via new non-convolution integral 
equations, which we are able so solve. See Theorems 4.2 and 4.3. 
In Section 5 we set up the integrals for the Cauchy problem for the 
more general equation Uu = Uxx + -|- Be~^)~^u. In Section 6 

we solve the integral equation in the case A = B = 1/2 and obtain 
a solution of the problem uu = Uxx + (sech^a;)M,x > 0, subject to 
u{x,D) = f{x),ut{x,0) = 0. (See Theorem 6.6). A brief conclusion 
follows. We believe that the methods of this paper may prove to be of 
use in the study of various types of linear PDEs. 

2. Introduction to the Method 

We illustrate the approach by solving the problem 

Uu = Uxx - 2 ^, u{x, 0) = f{x),Ut{x, 0) = g{x), x >0,t>0. 

oc 

This problem was solved using different symmetries in [S]. Here we 
use a simpler approach. We suppose that /, S' G >S(R+) the Schwartz 
space, consisting of smooth functions whose derivatives (including the 
function itself) decay at inhnity faster that any power, i.e. Schwartz 
functions are rapidly decreasing. It is known that the Fourier Trans¬ 
form and Mellin transforms are automorphisms of the Schwartz space, 
see [22]. This space is a topological vector space of functions (p such 
that G C'°°(M+) and x°‘ip^^\x) 0, x —)■ oo, a,/? G Mq. 

It is elementary that the scaling transformation u{x,t) —t u{Xx, Xt) 
is a symmetry. The idea is to introduce a new solution by setting 

poo 

U{x,t)= / (p{X)u{Xx, Xt)dX. (2.1) 

Jo 

With separation constant equal to one half, separation of variables 
leads to the solution ui{x,t) = \/xIu{,x)e^, where v = |\/4A -|- 1 and 
A(x) denotes the usual modihed Bessel function, (see 9.6.10 of m)- 
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We denote the Bessel function of the hrst kind by Ju{z), (9.1.10 of PQ). 
Applying the symmetry and replacing A with iX, we obtain a solution 

poo 

u{x,t) = / ip{X)'/XxJi^{Xx) cos{Xt)dX 

Jo 

poo 

+ / ip{X)\/~XxJ^{Xx) sm.{Xt)dX. 

Jo 

We absorbed the terms in i into ip and 'ijj. The requirement that 
u{x, 0) = f{x) leads to the equation 


f{x) = / (p{X)'/^Ji,{Xx)dX. 

Jo 

We immediately recognise that / must be the Hankel transform of (p. 
The inversion theorem for the Hankel transform, (see [ 23 ]) yields 


a(A) = / f{y)V^Ju{Xy)dy. 

Jo 

For the condition Ut{x, 0) = g{x) we hnd 

poo 

/ X'ip{X)'/)^Jiy{Xx)dX = g{x), 

Jo 

which gives 

1 /‘°° 


A 


( 2 , 2 ) 


(2.3) 


Hence 


poo poo 

u{x,t) = / / f{y)y/xyXJ^{Xx)J^{Xy)cos{Xt)dydX 

Jo Jo 

poo poo 

+ / g{y)y/3yyJ^{Xx)J^{Xy)sm{Xt)dydX, 

Jo Jo 

solves our initial value problem. Thus we have obtained a solution 
of the initial value problem from a separable solution and a scaling 
symmetry. This solution can be simplihed by making use of the result 


sin(ca:) {ax) {hx)dx 


= 0 , 


2\/a6 
cos(/i7r 
\fab 


Pn- 


/ 2 - 1/2 


IP' + aJ — 

2ab 


0 < c < b — a,0 < a < b 
, b — a<c<b + a,0<a<b 


TT 


52 flS _ ^2' 

■Qm-1/2 (- ^ - h 6 + a<c,0<a<6, 
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which is valid for /x > —1. (See formula 6.672.1 of [S]). Here Py and 
Qy are Legendre functions, d- We dehne 


p{x,y,t)= / J^{Xx)J^{Xy)sm{Xt)dX 

Jo 

via the preceding integral and the solution can be written 


(2.4) 


u{x, t) = 


f{y)^/JJy^Ju{Xx)J^{Xy) cos{Xt)dydX 


+ 


V^9{y)p{x,y,t)dy. 


(2.5) 


The absolute convergence of the second integral justihes our use of 
Fubini’s Theorem. 

The integral XJy{Xx)Jy{Xy) cos{Xt)dX is divergent, so we cannot 
reverse the order of integration in the hrst integral in fl2.5p . However, 
we can dehne a solution in the distributional sense by noting that for 
each n 

fu Q 

/ XJ^{Xx)J^{Xy) cos{Xt)dX = — / J^{Xx)J^{Xy) sm{Xt)dX. 

Jo ^ Jo 

Taking the limit as n —)■ cxo we dehne the distribution A on iS(M_|_) by 
setting 

/ cxD a poo 

j Ju{.Xx)Jy{Xy)su).{Xt)dXdy. (2.6) 

This leads to a solution in the distibutional sense given by 
d 


u{x, t) = 


dt 


+ / V^3(y)p[x,y,t)dy. 


(2.7) 


In general this procedure requires us to solve one or more integral 
equations if we are to satisfy the initial data. Often this equation turns 
out to be a familiar integral transform with a known inversion formula, 
as happened with our example. In this paper we will show how it can 
lead to entirely new solvable integral equations. 

This method hints at a deep connection between Lie symmetries and 
classical harmonic analysis and indeed this is the case. The relationship 
between Lie symmetry analysis and harmonic analysis and representa¬ 
tion theory has been developed over the years, beginning with 0. II 
and [6]. See also [U], [ID], [21], |2D] and [TB] . 

The basic idea of this work is that for linear PDFs, it is often possible 
to realise the Lie symmetries as global representations of the underlying 
Lie group in the following sense. One has a Lie group Q, a represen¬ 
tation {(T,V) of Q and a mapping A ■. V ^ H where iL is a solution 
space of the equation. So A maps f E V to some solution of the PDF. 
Then if the action of the symmetries on solutions u = Af, is denoted 
by p{g)u for g E Q, the following relationship holds. 
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{pi9)Af){x) = {Aa{g)f){x), g eG, (2.8) 

where f ^ V 

In fact this relationship has been established for many classes of 
PDEs and the previous references detail some of the work done in this 
area. 

A consequence of all this is that the operator 

Vi(x) = / ip(e)p(expevi)u(x)de, 

Jn 

is essentially a group theoretic Fourier transform, which inherits prop¬ 
erties of the representations equivalent to p. It is therefore not sur¬ 
prising that this technique so often leads to the problem of inverting a 
standard integral transform. In this paper we consider some problems 
where the integral equations do not appear in the existing literature. 

3. The Equation uu = u^x + f(x)u 

We will determine the forms of the potential function / which allow 
non trivial symmetries. For an in depth discussion of the computation 
of Lie symmetries, see the book by Olver [TB]. We look for a vector 
held V = ^(x,t)dx + T(x,t)dt + (j)(x,t,u)du which generates the Lie 
symmetries of the PDE. An application of Lie’s Theorem (Theorem 
2.31 pl04, m) leads to the dehning equations 

0^* = 0- + /(a;)0 + /'(a;)e«, (3.1) 

where 

G Gxx T (‘2‘Gxu ^xx)^x Xxx'^t T Guu'^x T (Gu 2.^x)'^xx ‘2'TxUxt 
0 Gtt T (‘2Gtu Xff)uf -f" Guu^t T (Gu ‘2 tG)u^i. 

From the dehning equations fl3.1l) . we hnd the conditions 

(j) OiU T /3, Pr, 2(j)xu ^xxi 

2‘(put Pti Prxj Gu 2,Tt (j)u ‘2‘^x 
4>tt + (Gu - 2Tt)f(x)u = (pxx + f(x)(p + f(x)u^. 

The function /3 will be an arbitrary solution of the PDE. Nothing fur¬ 
ther can be said about it. From the hrst and last of these equations we 
have 

att = axx +f(x)^ + 2Ttf(x). (3.2) 

It is also easy to see that ^xx = itt and Tu = Pex- Consequently 2(j)ut = 0, 
and this means a is independent of t. Similarly we see that a is also 
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independent oi x. It is thus a constant. From fl3.2p we see that 

2ex + = 0- (3-3) 

This means that t) = e~^ where c has yet to be hxed. 

With G{x) = ^, we now require 

c"{t)G{x) = G"{x)c{t). (3.4) 

This is only possible if G"jG = A a constant. We then have three 
possibilities for c : 

{C\t + C2 

Cl sin(ci;t) + C 2 cos{ut) (3.5) 

cie‘^* + C2e"^h 

Clearly the three possible choices for G are 

{ Ax + B 

(3.6) 

A cos ixx + B sin ixx 

where A, B are arbitrary constants. Correspondingly, the possible 
choices for / are 


fix) = 


1 

(Ax+B)^ 

1 

_ 1 _ 

(A cos u!x+B sinijx)^ ' 


(3.7) 


We can now write down a basis for the Lie algebra of symmetries for 
each choice of /. In each case, the Lie algebra is isomorphic to s [2 © r, 
where r is the Lie algebra generating the additive Lie group M. 

We are interested in the second case. The hrst case is obtainable 
from our illustrative example by a linear change of variables. If 


a basis for the Lie symmetry algebra is 


vi = e‘^\Ae‘^^ + Be-‘^^)d^ + - Be-‘^^)dt, 

V 2 = e-‘^\Ae‘^^ + Be-‘^^)d^ - e-^\Ae^^ - Be-^^)dt, 
V3 = dt, V4 = udu, = (3du. 


The vector field V 3 clearly induces the symmetry u{x, t) -A- u{x, t+e). 
Let us obtain the symmetry arising from the hrst vector held. We have 
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to solve 


dx 

de 


= x(0) = x, 


^ t(0) = t. 

CiC 


Adding these equations gives 


d{x + 1) 
de 




Now setting Z = x +1, we solve this first order DE and get 

^U){x+t) 


^Lu{x+t) _ 


1 - 2Aea;e‘^("^+d' 


Next we solve 


d{x — t) 
de 


= 2Be-^^^-^\ 


With z = X — i, 2:(0) = x — t we obtain the expression 

ga;(£-7) ^ ^ 2BeUJ. 


(3.8) 

(3.9) 


(3.10) 


(3.11) 


(3.12) 


Taking logs we have a pair of simultaneous equations for x and t. This 
gives us 




1 — 2 Aea ; e ‘^(*+0 


(3.13) 


and 


= 


^ixi(x+t) 


^(g2a;x ^ 25ea;e‘^("^+d)(l - 2Aea;e‘^("^+d)' 


From these we obtain expressions for x and t. From this we conclude 
that if M is a solution of 

1 


Hit Uxx T 


M, 


then so is 


u{x, t; e) =u ( - In 


In 


(Ae"^^ + Be-‘^^y ’ 

g 2 a;x ^ 2Beue‘^^^+^^ 

1 — 2 Aea ; e ‘^(*+0 

^uj{x+t) 

+ 2i?ewe^<^+*>)(l - 2A<^iV€^^^+^y) 


(3.14) 


, (3.15) 


The vector field V 2 can be exponentiated by similar means. This is 
left to the interested reader. 
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4. The Case of one zero constant 


To begin our analysis we will first consider the special case when 
T = 0, i? = 1. The case B arbitrary can be obtained from this by an 
obvious change of variables. This leads to the equation 

Utt = X e M. (4.1) 

For convenience we will make the change of variables y = e^. This leads 
to the equation 

Utt = y'^Uyy + yUy T |/ > 0. (4.2) 

We will solve this subject to the initial conditions u{y,0) = f{x), and 
Ut{y,0) = g{y) where f,g are suitable functions. It is straightforward 
to see that a stationary solution is Uo{y) = Jo{y). Application of the 
symmetry fld.lSp under the change of variables shows that 

Ui{y,t-,X) = Jo{^/y^ + 2Xye^), (4.3) 

is also a solution. Since t ^ —t is a symmetry, it is clear that 

U 2 {y, t; X) = Jo{\/y‘^ + 2Xye-^) (4.4) 


is again a solution. We now consider 


'“(».*)= 2 / + Jo(\/!FT2yy^) d\, 


(4.5) 


where ip has suitable decay to guarantee convergence of the integral. 
This is easily seen to satisfy the PDF fl4.2p and the conditions 

roo 

u{y,0)= / ip{X)Jo{^/y‘^ + 2Xy)dX (4.6) 

Jo 

and Ut{y, 0) = 0. Our hrst task is then to solve the integral equation 

V^{X)MVy'^ + ‘JXy)dX = f{y). (4.7) 

We proceed by taking the Laplace transform in y of both sides. We 
suppose that ip G Li(M+). Via the known inequality for the Bessel 
function i/x|W(x)| < C, x > 0, where C is an absolutely positive 
constant, we have from fl4.7l) the estimate 

^ (4.8) 

Thus 



f{y)e ‘"^dy 


< <^117^11x1(8+) 


'•°o g-sy 




dy 


V^C\\ip\\L,iR+) 
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Now 10.2.6, p59 of pS] gives 


JoWv'^ + 2Ai/)e 


g-A(v/ s2+i_s) 

\/s 2 + 1 


(4.9) 


If i^(s) = /g°° fiy)e ^^dy then an application of the Laplace transform 
and Fubini’s Theorem, (justified by fl4.8p ). gives 


(p(A)- 


3 —A(\/s^+l—5) 


-dX = F{s). 


\/s 2 + 1 

Letting $ denote the Laplace transform of <p, we have 

t(7?TT-s) _ 

The substitution s = sinh k reduces this to 

<h(e“^) = cosh fcF (sinh fc). 
Putting A; = — In ^ gives 

1 


^(0 = nU + 7M"U 7-e 




2 ve 


Thus 


(7(A) = 

We have thus established the following. 




(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 


Proposition 4.1. The equation Uu = y‘^Uyy + yUy + y‘^u with u{y,0) = 
f{y) and Ut{y, 0) = 0 has a solution 


u{y,t) 



(A) 


X 


JoiVv'^ + 2Ai/e*) + MVy'^ + 2Ai/e-*) dA, (4.15) 


provided the integral converges. Here F is the Laplace transform of f. 


Our next task is to determine an expression for the inverse Laplace 
transform. An elementary approach is to assume that / is analytic. 
That is 

CXD ^ 

/(a) = E )j/“(0)a". (4.16) 

In this case 

CXD ... 

F(.)=^/(»)( 0 )- 

71=0 


(4.17) 
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SO that 




Thus 


where 


ne yj {i-ey+^ 

(^(A) = 5^/(-)(0kn(A), 

(e^ + l)( 20 ' 


n=0 


gniX) = y 


-1 


(4.18) 


(4.19) 


(4.20) 


L 

which can be explicitly computed for any n. 

We can also obtain an explicit solution via a double integral as fol¬ 
lows. Making the substitution s = sinh(ln,^) = | ~ | j 5 ^ > 1 in 

the equation 


$(\/s2 + l-s) = \/s2 + lF(s). (4.21) 

we obtain 

I4 (1) = (4 + l) e -'>(^-' t ) f ( 2 y ) dy , { > 1. (4.22) 

To proceed, we use two operational relations for the Laplace trans¬ 
form, which may be found in [19]. Specihcally, if g{s) = (£h)(s), then 
via the uniqueness property for the Laplace transform 


1 fl 
-9 


= C 




which is (21) on page 171 of 


^ +il^fs-^)=7: 


Jo \2yXy) h{X)dX 
. Further 


(s), 


(4.23) 


'(9-A)‘''V(2vV5vy)ft(A).^ 


ry 

- {y - X)-^/^I.i {2yX{y - A)) h(A)^/A dX 


(s), 


(4.24) 


from formula (27), again pl71 of [T^]. Hence, returning to fl4.22p . we 
cancel the Laplace transform via the uniqueness property and take into 
account that Ii{z) = I-i{z) to obtain for i/ > 0 



Jo {2\/Xy^ ip{X)dX = 


ry Ii 

(J-2A)- 


( 2 yA(!/-A)) 

VX{y-X) 


f{2X)dX. 

(4.25) 


However the left-hand side of 04.251) represents a variant of the Han- 
kel transform of index zero, which admits the symmetric inversion for¬ 
mula for (f whose Mellin transform (p*{s) G Li(cr), cr = {s G C : 
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Re(s) = 1/2}. (See the details in | 25 ], Ch.2, Section 2.2, Example 2.2). 
We remark here that the Mellin transform is dehned by the integral 

poo 

/*('S) = / f{x)x^~^dx. (4.26) 

Jo 

Its inverse transform under certain conditions is given by the integral 

1 p'y-\-ioo 

/(^) = / ris)x-^ds (4.27) 

27Tt J ^—ioQ 

See [23], [25] for details of the Mellin transform. 

Thus, inverting the Hankel transform, we arrive at the unique and 
explicit solution of the integral equation fl4.7p . namely. 


(y9(A) = 


,2y 


Jo (‘^\/^ 


h (a/m(2i/ - m)) 


\/u{2y 


u 


{y - u)f{u) dudy, 

(4.28) 


A > 0. We note, that the interchange of the order of integration in 
fl4.28p is impossible because the corresponding integral with respect to 
y will be divergent. 

Now we observe that the inequality 


1 < r(z/+ 1) 



V 

Iu{x) < coshx < 


(4.29) 


see [15]) leads to 

l‘2y R {^^u{2y-u)^ 
0 \/u{2y - u) 


{y - u)f{u) du 


y^y 


r2y 


< -e 
- 2 


\f{u)\du. (4.30) 


Zemanian constructed a Frechet space i 3 ^(M+) characterized by three 
properties: Every 0 G i 3 ^(M+) is a rapidly decreasing smooth function 
on R+, with an expansion of the form 

(p{x) = (oq + Cl 2 x‘^ + • • • + CL 2 kX^^ + R 2 k{x)') , (4.31) 

where 

\im{x-^D)'^{x->^-J^(j){x)), (4.32) 

and the remainder satishes {x~^R 2 k{x) = o(l) as x —)■ O’*". (See 
Lemma 5.2.1 of [2H])- Zemanian proves that the Hankel transform of 
index y is an automorphism on i 3 ^(M_|_), [26]. We see then that if 

yJ^ !o^ \fi.'J)\du G i3o(K+), then (p G iDo(IR+)- 

We summarize our results by the following result. 


Theorem 4.2. Let f be a given continuous function on R+ such that 
f{y) = 0{y~^^‘^), 1 / —)■ oo. Let (p G Li(R+) and its Mellin transform 
(p*{s) G Li{a). Then ip given by formula if-2^ represents the unigue 
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solution of the integral equation O- Moreover ifye'^F{y) G i 3 o(I^+); 
where F{y) = \f{u)\du, and 


/(A) = 


^cxD n2y 


Jo ( 21 /^ 


'0 JO 


11 {^^z{2y - z)^ 
x/z(2y - z) 


{y - z)f{z)dzdy, 

(4.33) 


then 


= o / /(A) \jo{Vy‘^ + 2Ai/e*) + Jo{^/y‘^ + 2Xye-^) 


is a solution of 


'U'tt y ^yy F y^y F y 


d\, 

(4.34) 

(4.35) 


satisfying the initial conditions u{y,0) = f{y),Ut{y,0) = 0. 

To obtain a solution satisfying u{y, 0) = 0 and Ut{y, 0) = g{y) let 


I / p 

'<Jjiy,t) = - '0(A) JoiVy'^ F 2Xye^) - JoiVd^ + 2Ai/e-*) 


Then w{y, 0) = 0 and 


roo A 2 /Ji(v 2 /^ + 
wtiy,o) = - 0(A) 


dX. 


We therefore seek to solve 

Ji + 2\y^ 


A0(A)- 


^Jy^ F 2Xy 


Vy" + 2A|/ 


dX = -g{y) = g{y). 

y 


dX. 

(4.36) 

(4.37) 

(4.38) 


This can actually be reduced to the case covered by Theorem 14.21 In 
fact, assuming that A0(A) is absolutely continuous on M+ and its deriv¬ 
ative 93 (A) = [A0(A)]' satishes conditions of Theorem 14.21 we integrate 
by parts in fl4.38p . employing the equality 


y 


d 


and obtain the integral equation 

poo 

/ [A 0 (A)]'Jo(a/ 2/2 2Xy)dX = -g{y). 

Jo 

Hence our previous result gives 


0(A) = 


1‘QO I‘2y 

A Jo Jo Jo 


Jo ( 2 v^) 


h (^Vu{2y - u)^ 


u — y)g{u) dudydv, 
(4.39) 
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for A > 0. 

Theorem 4.3. Suppose that ye^G{y) G ioo(®+); G{y) = \g{z)\dz. 

Then the equation uu = y‘^Uyy+yUy+y‘^u with u{y, 0) = 0 and Ut{y, 0) = 
g{y) has a solution 


u{y,t) = - fif(A) Jo(V + 2Ai/e*) - Jo(\/|/^ + 2Ai/e *) 


dX, 

(4.40) 


where 


qW = 


1‘QO I‘2y 

^ Jo Jo Jo 


Jo i2y/yv) 


h [\/z{2y - z)) 
\/zi2y - z) 


z — y)g{z)dzdydv. 

(4.41) 


Hence, provided that / and g satisfy the conditions of Theorems 14.21 
and 14.31 a solntion of onr original problem may be written 

1 


u{y,t) = 


/(A) JoiVy'^ + 2Xye*) + Jo(\/y‘^ + 2Xye *) 


+ / aW + 2Xye^) - Jo{\/y'^ + 2Xye-^) 


dX 

dX ) . 

(4.42) 


4.0.1. The Direct Laplace Transform Approach. It is worth showing 
how the Laplace transform method can be applied to fl4.38p . This 
approach is compntationally nsefnl. On page 60 of [IH] we hnd 




(4.43) 


Jo \/y‘^ + 2Xy "" A 

If the Laplace transform of g{y)/y is G, then taking Laplace transform 
of both sides of fl4.38p . we obtain 


V'(A)(1 - e 


-A(Vs^ + l-s)\ _ 


)dX = G{s). 


(4.44) 


Which is the same as 


T(\/s2 + 1 - s) - T(0) = G(s). (4.45) 


We have denoted the Laplace transform of "0 by T. Pntting s = sinh k 
and fc = — In gives 


vl/(O-vI/(0) = G 



So that 


V^(A) = vI/(0)5(A) + 




(A). 


(4.46) 

(4.47) 


The term \k(0)(5(A) is in principle arbitrary, bnt making a choice for 
will £x it, and in any case, it actnally makes no contribntion to the 
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solution since the two terms involving Bessel functions cancel when 
evaluated at A = 0. 

Once more, if 


and so 


where 


= Z igS'”!'))!'" 

(4.48) 

1 i\ 

n=0 


-mphp . '‘PiP 

(4.49) 

oo 

V>(A) = >h(0)5(A) +J]9<”'(0)h„(A) 

(4.50) 

n=0 


~ H (1 _ ^2'^n+l ■ 

(4.51) 


Since the delta function terms make no contribution, the solution 
can be written 


w{y,t) = -j g{X) Jo{^/y^ + 2Xye^) - Jo{^/y^ + 2Xye-*) dX, 

(4.52) 

with g(A) = Er=o9*”’(0)*;n(5). 

Finding examples using the direct Laplace transform approach is not 
difficult. 

Example 4.1. We consider the solution in the case when /(y) = 0, 
9{y) = \y‘^iF 2 {\-, |, 2; -^) - \y. We hnd that 

r Q{y) . -1 .. _ 


f ^e-^ydy = - \ ='^(V^^Tl-s)-'^( 0 ). ( 4 . 53 ) 

Jo y 1 + s + vs^-Fi 

Putting s = sinh k gives 


-1 
1 + e 


- = xl/(e-^) - 11/(0) 


and ^ = e ^ leads to 


1 + e 1 + e 


1 = 4/(0-4/(0), 


and so 


^(e) = + ^(0) -1. 


(4.54) 

(4.55) 

(4.56) 


The inverse Laplace transform of both sides gives 
<^(A)=e-' + (4'(0)-l)5(A). 


(4.57) 
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The natural choice here is to take "^(A) = e ^ and then \l/(0) — 1 = 0. 
The solution our method gives for the initial value problem is then 


u{y,t) 



+ 2A?/e*) 


+ 2A?/e-*) 


d\. 

(4.58) 


This integral does not seem to be previously known. Although it is 
easy to evaluate numerically, it is worthwhile to give an evaluation. 

In order to calculate fl4.58p explicitly, we substitute the series expan¬ 
sion for the Bessel function (see 9.1.10 of jl] ) inside the integral and 
change the order of integration and summation via the absolute and 
uniform convergence. Thus we obtain 

poo 

ui{y,t)= / e~^Jo{^/y‘^ 2\ye^)dX 

Jo 

°° ( _ 1 \k rc>o 

Hence, the latter series becomes 


OO 


E 



CXD 


E 


m=0 


(2e*)"* 

{k — m)\y^ 
(4.60) 


Making a straightforward substitution in the inner sum and then chang¬ 
ing the order of summation, we get 


E 


k=0 


(_l)fc^2fc 

A^k\ 


E 

m=0 


(2e‘)™ 

{k — m)\y"^ 




1 

m\k\ 


EE - 


y_e 

2 


r 

’ 4 


m\{k + rn)\ 


m=0k=0 

The latter double series can be rewritten, recalling the dehnition of the 
Bessel function. Hence the value of the integral Ui{y,t) is 


<‘i(!/.«) = E(-l)'-TW(y). (4.61) 

fc =0 

The series fl4.6ip has no simple expression, but can be written in 
terms of the Lommel function of two variables [2], 

Uu{x,y) = '^(-1)^ l-j J 2 k+u{y), x,y > 0. (4.62) 

fc=o 

Thus appealing to relation (5.7.5.1) in |2], we write fl4.6ip in the form 
ui{y,t) = Uo{-eHy,y) + iUi{-eHy,y), 





















16 


MARK CRADDOCK AND SEMYON YAKUBOVICH 


where i = \/—l and, correspondingly, the valne of the integral fl4.58p 
as 

..r.. [Uo{-eHy,y) - Uo{-e-Hy,y)] + [Ui{-eHy,y) - Ui{-e-Hy,y)] 

'^[Ui — 2 

(4.63) 

The interested reader can prodnce many examples. For instance, 
the choice f{y) = 0,g{y) = yJo{y) leads to V’(A) = 2cos(A), for which 
T(0) = 0. 

4.0.2. The Case A ^ 0 and B = 0. . 

Now let ns briefly consider the equation 

Utt = Uxx + X G M. (4.64) 

Without loss of generality we take A = 1. Now set y = e~^. This again 
produces the equation 

Utt = y'^Uyy + yUy + y"^u, y > 0. (4.65) 

The analysis then proceeds exactly as in the previous case. 

5. The Case When A and B are both nonzero 
N ow consider the equation u''{x) + = 0. Suppose that 

u, AB > 0 and let a = Then with (3 = 

«(^) = (l + Ae"“A/?; 2/?; - 0 + f , 

is a stationary solution of the PDF Uu = u^x + Here 

2 F 1 is Gauss’ hypergeometric function, (Chapter 15 of [T]). Applying 
the symmetry we obtain the new, non stationary solution 


K{x,t;e) = 


1 + 

1 + 2Aeue‘^C+t) 




2i"i /9,/3;2A 


-(1 + 


1 + 2AeueA^+-t) 


We assume that f3 is not a negative integer here. We can immediately 
construct an integral operator mapping functions to solutions by setting 


u{x,t) = ^ <f{e)[K{x,f,e) + K{x,-t]e)]de (5.1) 


for suitable (p. This will satisfy Ut{y, 0) = 0 and 

A „2UJX ^ ^ 


m(x,0)= / (p(e) 


l + #e^ 


1 + 2Aea;e‘^^ 


2Fi /3,/3;2/3; 


-(1 + 

1 + 2Aea;e‘^^ 


de. 

(5.2) 
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If we take 


m(x, t) = 


^^{e) [K{x, t;e) — K{x, —t; e)]de, 


(5.3) 


(5.4) 


we obtain a solution satisfying u{x, 0) = 0 and 

/•oo g 

ut{x,0) = j 'il:{e)—K{x,t;e)\t=ode. 

If we require Ut{x, 0) = g{x), we have an integral equation for 

We will not attempt to solve these difficult integral equations here. 
Rather we look at the special case a; = 1 and A = B = 1/2. Then we 
obtain the solution 

/ 5 X. /2ee*+^ + e2^ - 1\ ^ /2ee*+^ + - 1\ 

u{x,t) = CX. [ --j -j , 

for the equation 

utt = Uxx + {sech^x)u. (5.5) 

In this = I (\/5 — l) . Suppose we set Ci = 0, C 2 = 1. Then 
form a solution of the PDE by setting 

. .4 r .4^ ^2ee‘+" + e2--T 
u{x, t) = J (p(e)g, ( -- 


de. 


we can 


( 6 . 6 ) 


Different domains of integration can also be considered. The operator 
dehned by fl5.6p maps functions (f to solutions of fl5.5p . We will show 
how this operator can be used to solve the initial value problem 

Utt = Uxx + (sech^a;)M (5.7) 

u{x,Q)=f{x), (5.8) 

Mi(a;,0)=0. (5.9) 

We suppose that / lies in some suitable function space, which will 
be discussed below. Now we remark that standard Laplace transform 
results and formula 23(3) on page 337 of [IS] gives the result 


\/|'e ® Tsinhi (e *^coshx) 

, V^sechx 


e = Q, 


/2ee‘+^ + - 1 

I a2x j_ T 


+ 1 


Thus formula fl5.6p can be rewritten 

(e-i^cosha:) 


pco poo 


w{x, t) = 


7’(e) 


\/^kecha; 


^^d^de 


® Tsmhx 2 / 1 (e *^coshx) 

^0 V ^S6cnx 

in which d) is the Laplace transform of ip. We obviously require suitable 
decay for p in order to justify the use of Fubini’s Theorem. Thus 
from the symmetry solution we have obtained an operator which maps 
functions $ from a suitable test space to solutions of fl5.5p . 
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If we let ^ ^ and use linearity we obtain the solutions 

cos(e“*.^ sinha;)e“^ {e~^^ coshx) 
wi{x,t) = - " 


W2{x,t) = 


y/^sechx 

sin(e“*^sinha:)e“5 (e“*^ coshx) 


\/^hechx 

From these we obtain the solution of fIS.Sp given by 


(5.11) 

(5.12) 


/ TT r°° cos(e *fsinha;)e *fcosha;) 

-21 

^ r°° sin(e“*fsinhx)e“ 2 J 1 (e“^fcosha;) 

^(0^- -7^=^ - -d^- 

Q \/^S6CilX 

Notice that time reversal is also a symmetry and so 

77 r°° cos(e*£sinha;)e^J,,,i (e*f coshx) 

u{x,t) = J- / $(0-^==ii-de 


\/^sechx 

^ sin(e*^sinhx)e^ Jj,,i (e*^coshx) 

'^'(0 - " - di , 

V'^sechx 

is again a solution of the PDF. We will therefore let 


m(x, t) = 


1 nr 


2 V 2 


<^>(0 


cos(e*^sinhx)e2 i (e*^ coshx) 


i/^sechx 


cos(e *,^sinhx)e 2J i(e coshx) 


-v/^sechx 

It is straightforward to see that 




utt = Uxx + (sech^x)M, 


i/^sechx 


and 


Mi(x,0) = 0. 


(5.13) 


nr r°° cos(f sinhx)J„,i (f coshx) 
u{x,0) = J- / «I>(0--de (5.14) 


(5.15) 


If we absorb the ^/^ into $, which is after all arbitrary, we see that 
if we can solve the integral equation 


/(x)\/sechx = y — J <I)(^) cos(^sinhx)J^ (^coshx) d^, (5.16) 

for 4), with fx = z^+ then we have a solution of our Cauchy problem. 
This is not a standard integral equation and does not seem to appear 
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anywhere in the existing literature. The solution is more involved than 
the previous integral equations we solved. Consequently, we will discuss 
its solution in the following section. 


6 . Solving the Integral Equation 


Let X G M_|_ and consider the following non-convolution integral equa¬ 
tions of the hrst kind 



ip{^) cos(^ sinha:) cosh.x)d^ = f{x) 


( 6 . 1 ) 



ip{^) sin(,^ sinha;) cosh.x)d^ = f{x), 


( 6 . 2 ) 


where f{x) is a given function and (p(^) is to be determined and 
Ju{z), H G C is the Bessel function of the hrst kind [T]. Though only 
the solution of the hrst equation is needed for our analysis, we present 
the solution of both. 

The key ingredients to solve integral equations fib.ip and fl6.2p are 
the following integrals (see [H], 6.699(2) and 6.699(1), p.723, with 
A = s - 1), 


X® ^ cos{bx)J^{cx)dx = 


2.-ir 


cr(i -F 


hLl£)'| 
2 ^ 


f s + v s — u 1 6^^ 
(6.3) 


where c > b, —Re(H) < Re(s) < 3/2, r( 2 ;) is the Euler gamma func¬ 
tion, see Chapter 6 of pQ. We also have 


X® ^ sin(6x) J,^(cx)dx = 


2"6r 

s + u + 1 s-l-1 


X 2F1 


u 


’ 2’ -2 


6 


. (6.4) 


As we see the left-hand sides of 06.311 and 06.411 represent the Mellin 
transforms of the kernels in equations 06.ip and 06.211 . Recall that for 
two functions /, g the generalized Mellin-Parseval equality holds (see 
[23] . [25] for the details) 


f{xy)g{y)dy = — 


^7+2CO 


f*{s)g*{l — s)x "^ds, X > 0. (6.5) 


' 7 — 200 


The Fourier cosine and sine transforms of the integrable function / are 
dehned by the formulas 


see [23] . 


(i^{^7/)(x) 


r°° (cos{xy) 
\l 71 Jo \sin(xi/) 


f{y)dy, 


( 6 . 6 ) 
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This equation does not have a unique solution. However if we impose 
certain extra conditions on (p, uniqueness can be obtained. We present 
the general case hrst. 

Theorem 6.1. LetRe(z/) > —l/2,(p(^) G S'(M+), f{x) G Li(M+), 
and G S'(M+), where g{x) = {x^ + / (sinh“^ x). 

Let also the Mellin transform ip* of ip satisfy the condition 


VJ*(s)|s|-"e(^/2-5)PI g ^ g [0,7r/2), 

where cr = {sGC, s=| + zr, rG M} and 

heir) = ip*{l /2 + zr)r(l/2 — ir + u) + p>*(l /2 — zr)r(l/2 + ir + u) 

(6.7) 

= O(T^). 

r ^ 0. Then the integral equation h6.1\) has the following solutions 


t(0 = 


P(s) 


Airi J^T{1 — s + u) 
1 


C"ds+ ^ ^ X 


vrWTT 


,(y+l/2)u 


V2 


+ V2 


V sin 


TT 


z/ + - ) - ) + cos(7rz//2) 


X 


L\/2e 


cos 


(e"r') ) - =in (e“r') (1 + ds 


Xi 


sm 




X / y '^cos(|/sinhM) / cos(i/sinhx)(coshx)^ ^/(^) dxdydu, 

Jo Jo 

( 6 . 8 ) 

where ^ > 0 and depending on an arbitrary function p{s), which is odd 
on a, i.e. 


p{s) = —p{l - s), s = 1/2 + ir. 


(6.9) 


Proof. Let a = {sgC, s = |+ir, tgM}. Since via integration by 
parts in the integral fl4.26p it is not difficult to verify that the Mellin 
transform p>*{s) of G 5'(M+) belongs to Li{a), we have by the inver¬ 
sion formula fl4.27p 

<d(0 = ^ (6-10) 

Substituting this integral into fl6.ip . we can proceed if we may change 
the order of integration. But since ( see 9.1.7 of ra) Jv{x) = O(x^), X — )■ 
0 and Re(z/) > —1/2, this will be possible if we show (see, for instance, 
in [25], Section 2.1) that there exists a positive absolute constant C 
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such that any hxed x > 0 and for almost all Ei,E 2 > 1 and all r G M 
f cos(^sinhx) Ji.(^coshx)(i^ < C. 

'El 

The latter inequality follows, employing the asymptotic behavior of the 
Bessel function at inhnity (see 9.2.1 of PP) 

Jy(x) ~ \ — cos (x -—) , X —)■ CXO 


TTX 


2 } 


and from the uniform convergence with respect to r of the following 
integrals 

I'E2 


'El 


TTZ/ 


^ COS (^(sinh x ± cosh x) =F uu ) 


( 6 . 11 ) 


This fact can be verihed with the use of the mean value theorems and 
integration by parts, and we omit the details. Hence with the use of 
the Boltz formula. 


2 Fi(a, b, c; z) = {1 - z) “ 2 ^^! ( a, c - 6; c; 


( 6 . 12 ) 


(see 15.3.4 of IH), in the right-hand side of fl6.3p . the equation (16.1 p 
takes the form 


^cosh'x r r(i^) 

= fix)- 


2 27ri 


El 


1 — S + V V + s 1 


2 ; ji-silih'Ul (i3*(s)(is 

(6.13) 

To proceed, we make the the substitution s = l/2-(-ir in the integral 
fl6.13p to obtain 


TT cosh^ X 
2 271 


1 /2—■ir+h' \ 

2 ) 


,Fi ( 1 / 2 ^. 1 . _ ^ 


2l/2+iTp I 1,/2+v+iT \ 


X (p*(l/2-|-ir)(ir = /(x). (6.14) 

Next we split the integral as -h and note the symmetry 

2 Ei{a,b;c; z) = 2 Eiib,a;c] z). The substitution r —r in the hrst 
integral and the duplication formula for the gamma function, 

1 


r{2z) = 


-,22- 




^/^T{z)T{z + l/2), 


(6.1.18 of PP) leads to the new form of (16.ip 

/•oo 2E1 


hcir)- 


( l,2+.+ir ^ l/2+.-ir . 1. _ ^j^^2 ^ 




( 3/2+y+iT'\ p f 3/2+1/—ir 

V 2 ) y 2 


dr = 2^/2+^ 


fix) 


cosh^(x) 

(6.15) 
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Moreover, the integral fl6.15p converges absolutely by virtue of the 
Stirling asymptotic formula for the gamma function ( see 6.1.37 of m) 
and the behavior at inhnity in r of the hypergeometric function (cf. 
Theorem 1.12) 




^2 —) = 0(1),T ^ oo. 


In the mean time, the latter hypergeometric function has the represen¬ 
tation (cf. [2], Vol. 2, relation (2.16.21.1), [2l], formula (1.101)) 


2 F 1 


1/2+ iy + iT 1/2+ iy — iT 1 ,9 . 

— -, — -; - smh^ x = 

2 2 2 


23/2—t 


X 


/ l/2+i/+iT\ p / 1/2+u—iT 

V 2 [ 2 


cos{ysmh.x)Kir {y)dy, (6.16) 


Re(z/) > —1/2, where Kir{x) is the modihed Bessel function of the 
third kind or the Macdonald function [1], and it is the kernel of the 
Kontorovich-Lebedev transform [23]. 

Substituting the right -hand side of 06.161) into the left-hand side of 
06.15p . we change the order of integration by the Fubini theorem, which 
is justihed via the absolute convergent of the iterated integrals. Hence 
after the use of the duplication formula for the gamma-function and 
the substitution u = sinh x, we arrive at the following integral equation 



r°^ hc{T) Kir{y)y’^ cos ( m |/) ^^^ 
0 |r (I + + iT) 1^ 


271 


f (sinh ^ u) 
(u^ + iy/2 ’ 


(6.17) 


since r(z)r(x) = r(z)r(z) = |r(z)|2. 

Further, by the dehnition of hdr), the Stirling asymptotic formula 
for the Gamma function, the uniform inequality for the Macdonald 
function m 


Kir{y)\ < e ^^/^iFo(//cos(5), 


5 e 



and our assumption y)*{s)\s\~'^e^'"F-&)\s\ ^ Li((t), one can apply the 
Fourier cosine transform to both sides of the equality fl6.17l) . Hence, 
returning to the original variable, we arrive at the equation 


0 \r[l + u + iT)\ 



cos(|/sinhx) 


y 


- 1 / 2 + 1 / 


(coshx)^ ‘'f{x)dx 


(6.18) 
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The left-hand side of fl6.18l) is related to the Kontorovich-Lebedev 
transform [21]. Recall that for snitable fnnctions g, we have 

2 poo POO 

g{x) = — / / Tsmh.{TrT)Kir{x)Kir{y)g{y)dydT. (6.19) 

X Jo Jo 

If we denote the right hand side of flb.lSp by (f){y) and introdnce 

7r^/ic(r) 


HJt) = 


2rsinh(7rr) |T (| -|- z/ -|- ir) | 


2 ’ 


then 


2 

4>{y) = ^ Tsmh{7rT)Hc{T)Kir{y)dT. 


TT^ 


( 6 . 20 ) 


( 6 . 21 ) 


Bnt the conditions of the theorem permit ns to invert the Kontorovich- 
Lebedev transform, so by fl6.19p we have 

' 0(l/) 


Heir) = 


y 


-Kir{y)dy. 


( 6 . 22 ) 


Thus we obtain 


^c(t) = ^rsinh(7rr) 


TT^ 


r ( - + iy + iT 




X 


cos(i/sinhx) 


(6.23) 


In the right-hand side of fl6.23p . we employ the Parseval equality for 
the Fourier cosine transform [23], which is allowed via conditions of 
the theorem. Then using relation (2.16.14.1) in [2], Vol. 2 and simple 
substitutions, the latter equality fl6.23p becomes 

2 POO 


^ ’ ttx/F V 2 / 


FI - + iy + iT 


cos{tu) 


X 


cos(i/sinhM) /■°° cos(i/sinhx) ^ 

^ —r— ^f(x) dxdydu. (6.24) 

(cosha;)^-! ^ ^ ^ ^ ^ 




Jo r JO 

Now, recalling the definition of hdr) in terms of (p*, let us write the 
equation fl6.24p with respect to the variable s = 1/2 -|-zr, dividing hrst 
both of its sides by F(l/2 -|- z/ -|- ir) and appealing to the reflection 
formula for the Gamma function. After making the change of variable 
t = e“, we have 

F(1-s-Fz/) S ^/2 T { l-s + u ) 

^ - s) = ^ 


F(s + z/) 


+ 


0F r F 


1/2-s 




cos(P(<-l)) 

Jo (cosha;)'^-! ^ ^ ^ ’ 

(6.25) 


y 


1 / 2 + 1 / 
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where s = 1/2 + ir, Re(z/) > —1/2. Further, denoting the iterated 
integral in the right-hand side of fl6.25p by 


F{s) = 


8\/2 


TT 


'1 


+ 


f+2 


COs(i|/(t - i)) 


'0 


,1/2+iy 


Jo (cosha; 
equation 06.251) becomes 

r(l — s + u 


°° cos(|/sinha:) , , , 

X / ^^ j[x)dxdydt, 


ip*{sy- 


r(s +1/) 


+ y?*(l-s) = 


F{s) r(i — s + u) 

r (r 


(6.26) 


(6.27) 


Note that F{s) = F{1 — s). 

Hence, 

\ r(s + i/) , ^ F(i - s) r(s + z/) 

^ ^—7 + 7^ (5) = 


g-1/2 \ p ( 1/2-s 


r(^ r, , 


or 


V7*(l -S) +V7*(s) 


r(i - s Fu) 

r(i — s + u) F{i — s) r(i — s + u) 


r(s + u) 


r (r (i^ 


(6.28) 


(6.29) 


Therefore, adding equations 06.27P and 06.29p . we hnd 

r(l-s + z/) [F{1 - s) + F{s)] T{l-s + u) 


2ip*{l - s) + 2if*{sy- 


r(s + z/) 




(6.30) 


Hence 


2^*{l-s) - 


F{s) r(l — s F y) 

r (r 


r(s Fi') = -r(i - s + z/) 


X 


2^*{s) 


F{1 — s)r(s -I- z/) 


s- 1/2 


1 / 2 -s 

2 


If we denote the right hand side of the latter algebraic equation by 
p(s), we see that it satishes the parity property 06.9p . Choice of either 
(p*(s) or (^*(1 — s) hxes the other. Hence hnal solutions of the equation 
06.ip can be expressed in terms of the inverse Mellin transform 04.27P 


7^(0 = —^ f ' dsF — [ 

r(l - s + 1^) 4m Ja r (r I 


(6.31) 
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where ^ > 0 and p{s) is an arbitrary fnnction, which is odd on the line 
a, i.e. satishes fl6.9p . The integrals converge via the relation fl6.25p and 
properties of (p* (s) that we have assnmed. Onr hnal goal is to calcnlate 
the second integral in the right -hand side of 06.311) . To do this, we 
will employ the generalized Parseval eqnality of type 06.5p (see [16], 
Chapter 7, Theorem 23), becanse the gamma-ratio behaves on a as 


r(l - s + u) 

r((s-i/2)/2)r((i/2-s)/2) 


0 (|j,|i+EeM)_ |Ims| ^ oo, 


and the inverse Mellin transform generally does not exist since the cor¬ 
responding integral 04.27P may diverge. So, recalling 06.45p . we observe 
that according to conditions of the theorem, the iterated integral 

poo poo 

/ cos{y{t — l/t)/2) / cos{y sinhx) {cosh f{x) dxdy 

Jo Jo 

is a fnnction of (t — l/t)/2 from the Schwartz space. Therefore, dne 
to the absolnte and nniform convergence of the corresponding integral 
06.45p . F(1 — s) is analytic in the right half -plane Re(s) > 1/2. Hence, 
retnrning to (16.3ip and nsing the dnplication formnla for the gamma 
fnnction r(l — s -)- z/) in the latter integral, we change the contonr a on 
the right-hand inhnite loop L+, encircling the right-hand simple poles 
of the obtained gamma fnnctions in the nnmerator. Precisely, with a 
simple snbstitntion we hnd 


1 [ F(1 - g) P(1 - g + u) 

AmJ^T{{s-l/2)/2)T{{l/2-s)/2f 
_ 2^ r F{1 - 2s) P((l + u)/2 - g)P(l + u/2 

“2^^/^ Jl+ P(s - 1/4)P(1/4- s) 


(6.32) 


Meanwhile, the integral 

_^ f r((i +- g)r(i + u/ 2 - s) 

2m P(s-1/4)P(1/4-s) 


(6.33) 


can be calcnlated, nsing the Slater theorem [2|, Vol. 3, and it can be 
expressed in terms of the hypergeometric fnnctions 1 F 2 . Conseqnently, 
we obtain 


1 

2'Ki 

x( 


P((l -|- l')/ 2 — s)P(l -|- Vj2 


P(s- l/4)P(l/4-s) 


I- 


.7 z/ 3 z/ 3 

1-^21 ^ ^ 

5 z/ 1 z/ 1 

1-^21 ^ ^ 


^z-^ds 


sin(f(z/+|)) 

^^l+ul2 



sin(f(z/ + I)) 


X 


(6.34) 
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However, the right-hand side of fl6.34p can be simplihed and written 
in terms of the elementary trigonometric fnnctions, appealing to rela¬ 
tion (7.14.1.1) in |2], Vol. 3 and keeping in mind particnlar cases of the 
modihed Bessel fnnction of the hrst kind. Thns we, hnally, derive 


-u 

20F 

— sin 


r((i +1^)/2 - s)r(i u/2 - ^ 


r(s- i/4)r(i/4-s) 


cos 


-|- u sin 


H—^sin 


3 X / TT 

- + U ] COS 


Z/ + 


1^ + 


V + 


+ 


cos(7rn/2) 




+ 5 cos 


(t) 


(6.36) 


Hence, combining with fl6.45p . fl6.3ip . fl6.32p . after straightforward 
manipulations we arrive at the hnal form fIb.Sp of solutions for the 
integral equation fib.ip . completing the proof of Theorem 1. 

□ 


Some simplihcation may be made to our solution. 


Corollary 6.2. When Re(i/) > —1/2, v ^ 2?7, -|- 1/2, n G No? hhe 
solutions 116. Si) take the form 


V?(0 = 


P(s) 


Airi — s + u) 


f ^ds -|- 


^(1+1^) 


sin I ^(z/ -I- - ) I X 


4ic+i)n 


VTWTT 


X 


i + ^/2Yl sinf^f.+ lV^Ucosf^^ 


^/2 


2 7 


L\/2e 


U V2 


cos ( —-^ sin ( — — sin ( — ( 1 H-^ 


V2j\ 


-sin 


^ sin(e“n 1-^ 


V 2 7 


2^ 


X 


/ (sinhu -I- sinhx)^ Isinhu — sinh^T —- dxdu, 

Jo V ^ ' ' J (coshx)^-i 

(6.36) 

for f>0. 


Proof. The proof will be completed if we compute the inner integral 
with respect to y in 06.81) after the corresponding interchange of the or¬ 
der of integration. In fact, let |Re(z/)| <1/2. Then since f{x) G 
Li(M+), the integral with respect to x in 06.81) converges absolutely 
and uniformly for y > 0. Hence the interchange is possible, where 
the relatively convergent integral by y can be calculated via relation 
(2.5.3.10) in |2], Vol. 1. This leads to 06.36p . Moreover, the obtained 
integral with respect to x in the right-hand side of 06.36P converges 
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absolutely and uniformly for u, Re(z/) > z/q > —1/2 under the condi¬ 
tion G Li(M+), representing an analytic function of u in the 

domain D = {u E C : Re(z/) > —1/2, u 2n + 1/2, n E Mq}. Hence, 
by analytic continuation equality fl6.36p holds for a\\ u E D. □ 

6.1. Conditions for Uniqueness. We now turn to the question of 
uniqueness of solutions for the integral equation. flG.ip . One way to 
obtain this is to assume the the Mellin transform (p* (s) satishes a par¬ 
ticular condition on a. 


Corollary 6.3. Suppose that f and tp satisfies the conditions of The¬ 
orem lh.il and the Mellin transform of ip has the additional property 
that 

</9*(s)r(l — s + v) = p*{l — s)r(s -|- z^), s E a (6.37) 

i.e. this product is even on a. Then the integral equation (a has a 
unique solution 




F{s)T{s + z/)r(l - s + 


2r 


1/2-s 

2 


s-1/2 

2 


(6.38) 


where F is given by /jd.fSl ). 

Remark 6.4. Equality fl6.37|) can be rewritten, for instance, in terms of 
the modihed Laplace transforms of the function p. Indeed, recalling 
the Mellin-Parseval identity fl6.5p . equality fl6.37p can be written in the 
following equivalent form 


X 


- 1 - 1 / 


exp{—t/x)p{t)Fdt = x'" j exp{—xt)p{t)Fdt, x > 0. 


Proof. Using the condition on the Mellin transform p*, we hnd from 
equation 06.301) that 

[F(l - s) + F{s)]T{s + z/)r(l -s + n) 


p*{s) = 


AT 


1/2-s 


s-1/2 


F{s)T{s -I- z/)r(l — s A- u) 


2T 


1/2-s 


s-1/2 


(6.39) 


since F{s) = F(1 — s) on a. Inverting the Mellin transform we end up 
with the unique solution in the form. 


i 


F{s)T{s + z/)r(l - s -f- 


2r 


1/2-s 

2 


s-1/2 

2 


(6.40) 


□ 


A natural question is when 'ip*{s) = (p*(s)r(l — s -|- z/) = — s) 

on cr? 
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Lemma 6.5. Let ip he the inverse Mellin transform of 
1p*{s) = 99*(s)r(l — S + u) = 1p*{l — s). 

If g{y) = is even, then ip*{s) = ip*{l — s) for s E cr. 

Proof. If s G (T, then s = I + ir. Clearly we require 

/ OO POO 

x~^P+"ip{^x)dx = / x~^^‘^~"ip{x)dx, 


(6.41) 


for every r G M. We convert this to a Fourier transform by setting 
X = e^. This gives 

/ OO POO 

P^ytP{ey)e^/^ydy = / e-^^yiP{ey)e^>^ydy. (6.42) 

-OO j —OO 

So the Fourier transform of gpy) = '^(e^)e^/^^ is even. This happens 
precisely when g is even. □ 

We are then able to write down a solution of our Cauchy problem. 

Theorem 6.6. LefRe(z/) > G 5'(M+), f{x) G Li(M+), 

and y~^~^/‘^{Fcg){y) G iS(M+), where g{x) = {x^ + f (sinh“^ x). 

Let also the Mellin transform (p* of p satisfy the conditions 

(/p*(s)|s|-"e("/2-<5)|s| ^ ^ g [0,7r/2), 


99*(s)F(1 — s + n) = p*{l — s)F(s + 10 ), s E a (6.43) 

where cr = {s G C, s = | + ir, r G M} and 

heir) = p*{l/2 + ir)F(l/2 — ir + u) + p*{l/2 — ir)F(l/2 + ir + z/) 

(6.44) 

= 0(r2), 
r ^ 0. Define 

r U-s, 1 1 r'“(F(c-j)) 


F(s)= C”“i + 


Then the problem 


has a solution 


cos(r/sinhO 


Utt = Uxx + (sech^a;)^, 
M(a;, 0 ) = f{x), ut{x,t)) = 0 , 


. , 1 r 

u{x,t) = — / 

TT* ./n . 


(6.45) 


K{t,x,f) 


F{s)T{s + z/)F(l -s + n) 


1/2-s \ p ( 5-1/2 


f Msdf, 


(6.46) 
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where 




cos(e*^ sinha;)e 2 cosh a:) 


+ 


■v/|sechx 

cos(e“*.^ sinhx)e“5 (e“*^ coshx) 
\/^sechx 


(6.47) 


Proof. We only need to establish convergence of the hnal integral and 
this follows from the conditions of the theorem. Specihcally, since the 
Fonrier cosine transform maps Schwartz fnnctions to Schwartz fnnc- 
tions, and the Mellin transform in s of a Schwartz fnnction, with 
Im(s) = t, is Schwartz in t, the fnnction F has rapid decay on the 
line a and hence its inverse Mellin transform has rapid decay. This 
guarantees the convergence of the hnal integral. □ 


Simplihcation of this result may be possible, but we will not discuss 
this question here. 


6.2. The Second Integral Equation. We hnish by turning now to 
the solution of the integral equation fl6.2p . This can be established in 
the same manner as was used for the proof of Theorem 16.11 with the 
use of the integral fl6.4p and is given by our next result. 

Theorem 6.7. Let Re(i/) > —^,95(0 ^ ‘5(K+), f{x) G Li(M+), 
andy~’^~^^‘^{Fsg){y) G iS(M+), where g{x) = / (sinh“^( x)). 

Let also the Mellin transform ip* of ip satisfy the condition 

(p*(s)|s|-^e("/2-5)PI g ^ g [0,71/2), (6.48) 

where a = {s G C, s = | + ir, r G M} and 


hsir) = p>*{l/2 + ir)r(3/2 — ir + u) + (p*(l/2 — ir)r(3/2 + ir + u) 
= 0 {t^), r^O. (6.49) 

Then the integral equation fh.gl) has the following solutions 


PiO = 


pis) 


C'ds + 


V2 ^-(2+^^) 


X 


Am r(2 — s + n) 

\ 


TTx/Tr 


' —oo 
/*oo 


g(i^+3/2)« coshu 


X / yi/2 '^cos(i/sinhM) / sin(i/sinha;)(cosha;)^ ’^fi^) dxdydu, 

Jo Jo 

(6.50) 


for ^ > 0, depending on an arbitrary function pis), which satisfies 
relation on a. 
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Proof. In fact, an analog of equation fl6.13p will be the equality 
^sinh X cosh'^ x 


27ri 


2r(^) 


2-^1 


-—ds = f{x). (6.51) 


Hence, by a similar calculation to that leading to fl6.15p . we get 

poo 

/ hsir) 


2 F 1 (3/2^^ 3/2^. |. _ ^ 


f b/2+u+iT\ p f 5/2+1^—ir 

V 2 ) i 2 


dr = }(r) 


sinha; 


(6.52) 


and since 
sinha; 2-^1 


3/2 + i/ + ir 3/2 + 1 / —ir 3 . 2 


2 

^1/2—u 


r((3/2 + z/ + ir)/2)r((3/2 + v- iT)/2) 


; -sinh x ) = 


y"' ^/‘^sixi{ysmh.x)Kir{y)dy, 
(6.53) 


Re(z/) > —1/2, we derive an analog of the equation fl6.17l) 


'0 ao 


r(| + 1 / + ir)r(| + z/ - ir) 


(6.54) 


Further, according to conditions of the theorem we may take the 
Fourier sine transform and then the Kontorovich-Lebedev transform, 
to the function hs(r) in the form 

2 noo 


hs{T) = —rsinh(7rr) 


TT^ 


+ V + ir 




X / sin(?/sinhx)(cosha;)^ dxdy. 


(6.55) 


Moreover, after the use of the Parseval equality for the Fourier sine 
transform [23], relation (2.16.14.1) in [2], Vol. 2, integration by parts 
and differentiation under the integral sign in the integral with respect 
to M, which is permitted under conditions of the theorem, we obtain 

2 poo poo poo 




TT+TT 


F 


+ u + ir 


cos{tu) coshu 


'0 JO JO 


X y^^'^ ^ cos(j/sinh u) sin(?/sinh x) (cosh a;)^ f{x) dxdydu. 


(6.56) 


Then similarly to equations fl6.25p . fl6.45p and fl6.27p we derive from 
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^ ^+ ^ (1 - s) = 


G'(s) r( 2 -s + z/) 


r(i + s + 1 /) 


r((s + l/2)/2)r((3/2-5)/2) 


, s e cr, 
(6.57) 


where 


G(.) = 


2^/2 


r*oo /*oo /*oo 


TT 


g“(i/2 *) cosliM ^cos(i/sinh 




' —OO ^0 «/ 0 


X sin(i/sinha;)(cosha;)^ dxdydu, 


(6.58) 

Hence as in the proof of Theorem 16.11 the hnal solutions of the equation 
fl6.2p can be expressed in terms of the inverse Mellin transform fl4.27p 




G'(l-s) r(2-s + z/)t 


^+1 

2 


I- 


-ds. 


(6.59) 


Then as above, 

If G(1 - s) r(2 - s + p) 
4^Xr((s + l/2)/2)r((3/2-s)/2) 


r^ds = 


2^+1 f ^(1 - 2s) r(l + Z//2 - s)r(3/2 + Z//2 - s) 


27riy/7r 


r(s + l/4)r(3/4-s) 


i^erds. 

(6.60) 


Meanwhile, the integral 
1 


" r(l + I//2 - s)r(3/2 + I//2 - s) _ 

2m r(s + l/4)r(3/4 - s) ^ ^ 

can be calculated in terms of the elementary functions. Precisely, we 
derive 


r(i + |-s)r(| + f-s) _ 


sm 


27ii 


r(s + i)r(f-s) 


2 ; ^ds = 


(^+ 7 ) - 




1 +- 
2 


(6.61) 

Hence, combining with fl6.58p . fl6.59p and fl6.60p . we arrive at the hnal 
form fl6.50p of solutions for the integral equation fl6.2l) . completing the 
proof of Theorem 16.71 

□ 

The corresponding corollary for the values z/gC: Re(z/) > —1/2, u ^ 
2n + 7/2, n G No can be formulated as follows 

Corollary 6.8. When Re{p) > —u 7 ^ 2n+7/2, n G Nq, the solutions 
H6.50\) take the form 


AO = 


pO) 


Am r(2 — s -\r p) 


f-’^ds+f-^‘^+'^A{?>/2-p) sin( 7 r( 3 / 2 -z/)/ 2 ) 
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X 


7r\/^ J- 

OO 


X 


(^+3/2)«gj^ cosh 

(sinhw + sinhx)^”^^^ (coshx)^“‘^/(x) dxdu 


u 


uo 


— / (sinhu — sinhx)*^ ^^^(coshx)^ '^f{x)dxdu 


+ / (sinhx — sinhw)*^ ^'^^(coshx)^ dxdu 


, e>o. 


u 

We conclude with a condition that will guarantee unique solutions 
for the second integral equation fl6.2lh The proof is similar to that for 
Corollary 16.31 and we omit it. 


Corollary 6.9. Let f and ip satisfy the conditions of Theorem \6. 7| and 
suppose further that 


</9*(s)r(2 — s + u) = (y9*(l — s)r(l + s + u), s G a. 
Then equation 46'. iH) has a unique solution given by 
1 


= 


471* Ja p 3/2 -s I Y 

where G is given by ^6f53^. 


G(i - i,)r(i + s + ^)r(2 - s +_0 (g g2j 




7. Conclusion 

The methods introduced in [8] and extended in the current work lead 
to many interesting problems involving integral transforms and integral 
equations. There are still very many open questions, even in terms of 
the focus of this paper. For example, one would like to be able to solve 
the general integral equation fl5.2p . When applied to other PDFs, the 
method also generates integral equations that do not seem to have been 
studied in the literature and many interesting questions arise from the 
investigation of these problems. We hope that this work stimulates 
further research in this area. 
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